A new transiting extrasolar giant planet 

Maciej Konacki*, Guillermo Torres^, Saurabh Jha'* & Dimitar D. 
Sasselov^ 

* Div. of Geological & Planetary Sciences 150-21, California Institute of Technology, Pasadena, CA 91125, 
USA 

t Harvard- Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 

* Department of Astronomy, University of California, Berkeley, CA 94720, USA 



A conceptually simple and technologically feasible approach to finding plan- 
ets orbiting other stars is to observe the periodic dimming of starlight due 
to a planet transiting in front of its star. Despite many intense photomet- 
ric searches, no transiting planet had yet been discovered in this way. The 
only known transiting extrasolar planet, 1,2 HD 209458b, was first detected by 
precise radial velocity measurements. 2 ' 3 We have measured radial velocities 
of a star, OGLE-TR-56, which shows a 1.2-day transit-like light curve found 
photometrically by Udalski et al. 4,5 Here we show that the velocity changes we 
detect are probably induced by an object of 0.9 Jupiter masses — a very close- 
in gas-giant planet only 0.023 AU from its star, with a planetary radius of 1.3 
Jupiter radii and a mean density of ~0.5 g cm 3 . At its small orbital distance, 
OGLE-TR-56b is hotter than any known planet, approaching 1900 K, but 
it is stable against long-term evaporation or tidal disruption. As the planet 
with the tightest known orbit, OGLE-TR-56b will place strong constraints on 
planet formation and migration models. 



The advent of high-precision Doppler and timing techniques in the past decade has 
brought a rich bounty of giant planets 6-8 as well as smaller, terrestrial-mass pulsar plan- 
ets. 9 To date over one hundred extrasolar giant planets have been found by different 
groups using precise radial velocity measurements. 8 Photometric observations of tran- 
siting planets, when combined with radial velocities, yield entirely new diagnostics: the 
planet size and mean density. 1 ' 2 Transits supply the orbital inclination and a precise mass 
for the planet, and they additionally enable a number of follow-up studies. 10-13 Hence, a 
large number of transit searches are already underway or under development. 14 However, 
photometry alone cannot distinguish whether the occulting object is a gas giant planet 
(~1-13 Jupiter masses), a brown dwarf (~13-80 Jupiter masses) or a very late type dwarf 
star, because such objects have nearly constant radius over a range from ~0.001 to 0.1 
Solar masses. This critical parameter, the mass of the companion, can be determined 
from the amplitude of the radial velocity variation induced in the star. 

One of the most successful searches to date is the Optical Gravitational Lensing Exper- 
iment (OGLE), which uncovered 59 transiting candidates in three fields in the direction 
of the Galactic centre (OGLE-III), 4 ' 5 with estimated sizes for the possible companions of 
~l-4 Jupiter radii. The large number of relatively faint (V = 14-18 mag) candidates 
to study led to our strategy of a preliminary spectroscopic reconnaissance to detect and 
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reject large-amplitude (high-mass) companions, followed by more precise observations of 
the very best candidates that remained. Of the 59 OGLE candidates, 20 were unsuitable: 
one is a duplicate entry, 4 have no ephemeris (only one transit was recorded), 8 show ob- 
vious signs in the light curve of a secondary eclipse and/or out-of-eclipse variations (clear 
indications of a stellar companion), and 7 were considered too faint to follow up. We un- 
dertook low-resolution spectroscopy of the other 39 candidates in late June and mid- July 
2002 on the Tillinghast 1.5-m telescope at the F. L. Whipple Observatory (Arizona) and 
the 6.5-m Magellan I Baade telescope at Las Campanas Observatory (Chile). These spec- 
tra were used to eliminate stellar binaries, which can produce shallow, planet-like eclipses 
due to blending with light from another star, grazing geometry, or the combination of 
a large (early-type) primary and a small stellar secondary, but are betrayed by large, 
easily-detected velocity variations (tens of km s" 1 ). We found 25 of the 39 candidates 
to be stellar binaries, and 8 to be of early spectral type. Only 6 solar-type candidates 
remained with no detected variations at the few km s _1 level (G. Torres et al., in prep.). 

Subsequently, we used the high resolution echelle spectrograph (HIRES) 15 on the 
Keck I 10-m telescope at the W. M. Keck Observatory (Hawaii) on the nights of July 
24-27, 2002, to obtain spectra of 5 of these candidates and measure more precise veloci- 
ties. OGLE-TR-3 turned out to be the result of grazing eclipses and blending (with even 
a hint of a secondary eclipse present in the light curve), and the data for OGLE-TR-33, 
OGLE-TR-10, and OGLE-TR-58 are as yet inconclusive and require further measure- 
ments. OGLE-TR-33 exhibits a complex spectral line profile behaviour and could also 
be a blend. OGLE-TR-10 shows insignificant velocity variation, which is consistent with 
a sub- Jovian mass planetary companion; OGLE-TR-58 is still inconclusive because of 
the uncertain ephemeris (M. Konacki et al., in prep.). Only OGLE-TR-56 showed clear 
low-amplitude velocity changes consistent with its 1.21190-day photometric variation, 5 
revealing the planetary nature of the companion. With only one bona-fide planet (or at 
most 3) among the 39 + 8 objects examined spectroscopically or ruled out on the basis of 
their light curves, the yield of planets in this particular photometric search has turned out 
to be very low: at least 94% (possibly up to 98%) of the candidates are "false positives". 
This is likely to be due in part to the crowded field towards the Galactic centre, which 
increases the incidence of blends. 

We report here our results for OGLE-TR-56 (I ~ 15.3 mag). Radial velocities were 
obtained using exposures of a Th-Ar lamp before and after the stellar exposure for wave- 
length calibration, and standard cross-correlation against a carefully-matched synthetic 
template spectrum (see Table 1). Our nightly-averaged measurements rule out a constant 
velocity at the 99.3% confidence level, and are much better represented by a Keplerian 
model of an orbiting planet (Figure la and lb). Note that the period and phase of the 
solid curve are entirely fixed by the transit photometry, as the ephemeris is constrained 
extremely well by the 12 transits detected so far (A. Udalski, private communication). 
The only remaining free parameters are the amplitude of the orbital motion (the key to 
establishing the mass of the companion) and the centre-of-mass velocity, both of which 
can be accurately determined from our velocity measurements. The properties of the 
planet and the star are summarised in Table 2, and Figure lc shows a phased light curve 
of the transit together with our fitted model. 

We performed numerous tests to place limits on any systematic errors in our radial 
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velocities and to examine other possible causes for the variation. These are crucial to 
assess the reality of our detection. On each night we observed two "standards" (HD 209458 
and HD 179949) which harbor close-in planets with known orbits. 3 ' 16 We derived radial 
velocities using the same Th-Ar method as for OGLE-TR-56, and also using the I2 gas 
absorption cell to achieve higher accuracy than is possible for our faint OGLE candidates. 
In Figure 2 we show that the measured velocity difference between our two standards 
(HD 209458 minus HD 179949) is similar for the Th-Ar and I 2 techniques, and more 
importantly, that both are consistent with the expected velocity change. This indicates 
that we are able to detect real variations at a level similar to those we see in OGLE-TR-56. 

We can rule out the possibility that OGLE-TR-56 is a giant star eclipsed by a smaller 
main-sequence star, both from a test based on the star's density inferred directly from 
the transit light curve, 17 and from the very short orbital period. We also examined the 
spectra for sky/solar spectrum contamination from scattered moonlight; a very small 
contribution from this source was removed using TODCOR, a two-dimensional cross- 
correlation technique. 18 The separation between the sky lines and the stellar lines is large 
enough (~30 km s" 1 ) that the effect on our derived velocities is very small. 

Blending of the light with other stars is the most serious concern 19,20 in a crowded 
field such as toward the Galactic centre. We have examined the profiles of the stellar 
spectral lines for asymmetries and any phase-dependent variations that can result from 
blending. Very little asymmetry is present, and no correlation with phase is observed. 
In addition, we performed numerical simulations to fit the observed light curve assuming 
OGLE-TR-56 is blended with a fainter eclipsing binary. Extensive tests show that with 
a photometric precision similar to the OGLE data (er ~ 0.003-0.015 mag), almost any 
transit-like light curve can be reproduced as a blend, and only with spectroscopy can 
these cases be recognized. For each trial simulation, the relative brightness and velocity 
amplitude of the primary in the eclipsing binary can be predicted. Although a good fit 
to the photometry of OGLE-TR-56 can indeed be obtained for a model with a single 
star blended with a fainter system comprising a G star eclipsed by a late M star, the G 
star would be bright enough that it would introduce strong line asymmetries (which are 
not seen), or would be detected directly by the presence of a second set of lines in the 
spectrum. Careful inspection using TODCOR 18 rules this out as well. Therefore, based 
on the data available, a blend scenario seems extremely unlikely. 

This is the faintest (V ~ 16.6 mag) and most distant (~ 1500 pc) star around which 
a planet with a known orbit has been discovered. The planet is quite similar to the only 
other extrasolar giant planet with a known radius, HD 209458b, except for having an 
orbit which is almost two times smaller. Thus its substellar hemisphere can heat up to 
~1900 K. However, this is still insufficient to cause appreciable planet evaporation (with 
thermal r.m.s. velocity for hydrogen of ~7 km s _1 compared to a surface escape velocity of 
~50 km s _1 ). The tidal Roche lobe radius of OGLE-TR-56b at its distance from the star 
is ~2 planet radii. The planet's orbit is most likely circularised (e = 0.0) and its rotation 
tidally locked, but the star's rotation is not synchronised (vsini ~ 3 km s _1 ). Thus the 
system appears to have adequate long-term stability. Interestingly, OGLE-TR-56b is the 
first planet found in an orbit much shorter than the current cutoff of close-in giant planets 
at 3-4 day periods (~0.04 AU). 8 This might indicate a different mechanism for halting 
migration in a protoplanetary disk. For example, OGLE-TR-56b may be representative of 
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a very small population of objects — the so-called class II planets, which have lost some 
of their mass through Roche lobe overflow 21 but survived in close proximity to the star; 
a detailed theoretical study of OGLE-TR-56b will be presented elsewhere (D. Sasselov, 
in prep.). These observations clearly show that transit searches provide a useful tool in 
adding to the amazing diversity of extrasolar planets being discovered. 
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Date 


RV 


Error 


MJD 


kins" 1 


kins- 1 


52480.4239 


-49.26 


0.20 


52481.4011 


-49.44 


0.08 


52481.4178 


-49.24 


0.09 


52483.3984 


-49.60 


0.06 


52483.4152 


-49.78 


0.11 



Table 1. OGLE-TR-56 radial velocities. The velocities (reduced to the solar system 
barycentre) and formal errors are given for each of our individual spectra of OGLE-TR-56. 
The data indicate a significant variation; a flat line fit gives x 2 — 20 with 4 degrees of 
freedom (0.06% false alarm probability), which is considerably worse than a fit to a 
Keplerian orbit model with a fixed ephemeris (x 2 — 5 with 3 degrees of freedom; 17% 
probability). Having shown, as a check, that the velocities from separate exposures on 
the same night (originally intended for cosmic ray removal) are not significantly different, 
we have adopted the nightly averages for subsequent use. A similarly high significance is 
found for the conclusion that the average velocities are not well fit by a flat line (99.3% 
confidence level). 



Parameter 


Value 


Velocity amplitude 


0.167 ±0.027 km s" 1 


Centre-of-mass velocity 


-49.49 ±0.02 km s" 1 


Orbital period 


1.21190 ±0.00001 days 


Reference transit epoch (MJD) 


52072.185 ±0.003 


Star mass 


1.04 ±0.05 M 


Star radius 


1.10 ±0.10 R 


Limb darkening coefficient (/ band) 


0.56 ±0.06 


Orbital inclination 


86 ± 2 deg 


Planet distance from star 


0.0225 AU 


Planet mass 


0.9 ± 0.3 M Jup 


Planet radius 


1.30 ±0.15 R Jup 


Planet density 


0.5 ± 0.3 g cm" 3 



Table 2. Derived stellar and planetary parameters. The physical properties of 
the star were derived by modelling the high-resolution spectra with numerical model 
atmospheres. We find that OGLE-TR-56 is very similar to the Sun, with a temperature 
of T eS ~ 5900 K. The star's mass and radius were computed from our stellar evolution 
tracks. 23 Combining the stellar parameters with the OGLE-III /-band photometry yields 
the planetary radius and orbital inclination. The uncertainties shown for the orbital 
elements are formal errors; the errors for the planet mass and radius additionally reflect 
our conservative estimate of systematic uncertainties. 
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Figure 1. Spectroscopic and photometric observations, (a) Our radial velocity 
measurements for OGLE-TR-56 (nightly averages). Note the good agreement with a 
curve whose only free parameters are the amplitude and systemic velocity; the period and 
transit epoch are fixed from the OGLE-III photometry. The curve assumes a circularised 
orbit (e = 0.0), as is theoretically expected, and the line thickness corresponds to the 
phase uncertainty, (b) % 2 surface showing the confidence region for our determination of 
the velocity amplitude and centre-of-mass velocity. The detection of a change in velocity 
(nonzero velocity amplitude) is formally significant at 6a (see also Table 2). (c) The 
photometric transit light curve of OGLE-TR-56 from Udalski et al. 5 The transit has 
an extended flat bottom, and its 1.2%-depth points to a Jupiter-size body (given our 
determination of a Sun-like primary). The solid line represents our fitted solution to 
derive the system parameters. The light curve shows no evidence of other variations or 
of a secondary eclipse (the hallmark of a strongly blended stellar binary). 
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Figure 2. Tests for systematic errors. Predicted radial velocity difference between 
our two standard stars with known planets — HD 209458 3 and HD 179949 16 (solid line), 
compared with measurements on each of our four observing nights in July 2002. The 
filled circles are our Th-Ar velocity differences (HD 209458 minus HD 179949, from the 
blue echelle orders beyond the iodine spectrum cutoff) with a typical internal uncertainty 
of about 100 m s _1 . These differences should be independent of the wavelength solution 
itself, and should reveal only the real difference in the Doppler shifts of the stars as well 
as any systematic problems of an instrumental nature. For comparison, our more precise 
iodine-cell velocity differences for the same stars (squares) have typical uncertainties of 
20 m s~ 1 . The uncertainty introduced by errors in the orbital elements of the standards 
is indicated by the shaded area. The graph shows that we have succeeded in measuring 
small changes in velocity on different nights using the standard Th-Ar technique, which 
reflects on the excellent stability of the HIRES instrument. The same technique was 
applied to our observations of OGLE-TR-56. 



